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Abstract. We consider a sequence of polynomials appearing in expressions 
for the derivatives of the Lambert W function. The coefficients of each polyno- 
mial are shown to form a positive sequence that is log-concave and unimodal. 
This property implies that the positive real branch of the Lambert W function 
is a Bernstein function. 



1. Introduction 

The Lambert W function was defined and studied in [5] • It is a multivalued func- 
tion having branches Wk{z), each of which obeys Wkexp(Wk) = z. The principal 
branch Wn maps the set of positive reals to itself, and is the only branch considered 
here. Therefore we omit the subscript for brevity. The nth derivative of W is 
given implicitly by 

d n W(x) eM-nW(x)) Pn (W(x)) 

(1-1) — ; = t, TTr/ », — ; tor n > 1 , 

v ; dx n (1 + W(x)) 2 ' 1 - 1 

where the polynomials p n (ui) satisfy pi(w) — 1, and the recurrence relation 
(1.2) p n+ i(w) — —(nw + in - l)p n (w) + (1 + w)p' n (w) for n > 1 . 
In [B], the first 5 polynomials were printed explicitly: 

pi(w) = 1 , P2(w) = —2 — w , P3(w) = 9 + 8w + 2w 2 , 

p 4 ( w ) = -64 - 79w - 36w 2 - 6w 3 , 

p 5 (w) = 625 + 974w + 622w 2 + 192w 3 + 24w 4 . 

These initial cases suggest the conjecture that each polynomial (— \) n ~ 1 p n (w) has 
all positive coefficients, and if this is true, then dW(x)/dx is a completely mono- 
tonic function [11] , We here prove the conjecture and prove in addition that the 
coefficients are unimodal and log-concave. 

2. Formulae for the coefficients 
In view of the conjecture, we write 

n-l 

(2.1) p n (w) = (-l) n - X J2Pn,kW k ■ 

k=0 

We now give several theorems regarding the coefficients. 
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Theorem 2.1. The coefficients (3 n ,k defined in (12. lj) obey the recurrence relations 

(2.2) /3„, = n"" 1 , p nA = 3n n - (n + 1)" - n n ' Y , 

(2.3) /?„,„_! = (n - 1)! , j 9 n , n _ 2 = (2n-2)(n-l)! > 

(2.4) /3 n+ i, fc = (3n-A;-l)/3 n , fc + n/3 n , fc _i-(A; + l) j 8 ni fc + i , 2 < fc < ?i - 3 . 
Proof. By substituting (|2.1I) into (|1.2p and equating coefficients. □ 

Theorem 2.2. An explicit expression for the coefficients /3 n ,k is 

m=0 v 7 q=0 v y 7 



m+n— 1 



Proof. We rewrite (jl.ll) in the form 

Pn(W(aO) = (1 + ^(x)) 2 "-^"^^)^^ . 

From the Taylor series of W(x) around x — 0, given in [5], we obtain 

d n W(x) = y ( z mT^ a;m _„ _ 
da;™ ' (m — n)! 

m—n 

Substituting this into the expression of p n , using x = We w and changing the index 
of summation, we obtain the equation 

OO Q 

(2.6) pjw) = (1 + w) 2 "" 1 V(-l) n+s - x (n + s )»+«~i^ e («+^)- . 

We expand the right side around w — and equate cocfhcicnts of w. □ 



Remark 2.3. The polynomials p n (w) can be expressed in terms of the diagonal 
Poisson transform D n [/ S ; z] defined in [TU], namely, by (|2.6[) 

(2.7) p„M = (-l)"-^! + w) 2( "- 1) D n [(n + a)"" 1 ; -«;] . 

Theorem 2.4. 77ie coefficients can equivalently be expressed either in terms of 
shifted r -Stirling numbers of the second kind {^_^.} defined in [3], 

(2-8) ^ = E(-ir( 2 ; _1 )( 2n " 1+w } , 

m=0 v 7 k J n 

or in terms of Bernoulli polynomials of higher order Bn \\) defined in 



m=0 



or m terms of the forward difference operator A [71 p. 188], 



,fe — 

m=0 



''' 7 2n-l\(-l) 



• / -^ ' fc — to / m! 
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Proof. We convert (|2.5|) using identities found in [3] and [8] respectively. 

i (q + r) n 



and 



(Ti = i rE(- i r- 

I m + r I „ to! ^— ' 

(m + raj! ^— ' \q/ 



g=0 

□ 

3. Properties of the coefficients 

We now give theorems regarding the properties of the f3 n ,k- We recall the follow- 
ing definitions [12] . A sequence Co, c\ . . . c n of real numbers is said to be unimodal 
if for some < j < n we have Co < c\ < . . . < c,- > c J+ i > . . . > c„, and it is 
said to be logarithmically concave (or log-concave for short) if c^-iCk+i < c? for 
all 1 < fc < n — 1. We prove that for each fixed n, the /3 n> fc are unimodal and 
log-concave with respect to k. Since a log-concave sequence of positive terms is 
unimodal 15J, it is convenient to start with the log-concavity property. 

Theorem 3.1. For fixed n>3 the sequence {k\f3 n ,k}2 = o * s log-concave. 
Proof. Using (|2.5p we can write 



kW n ,k = (2n- 1)! ( k ^x m yk-m 

m—0 



where 



(3.1) x m = y( m )a J , aj = (-iy(n + j) m+n -' 



3=0 



and y m — l/(2n — 1 — to)! . Since the binomial convolution preserves the log- 
concavity property [13 ] 114 ] . it is sufficient to show that the sequences {cc m } and 
{y m } are log-concave. We have 

aj-ia j+l = (-ly-^n + j - l) m +"- 1 (-l)-'"+ 1 (n + j + 

= (-l) 2 ^ ((71 + jf - I)"*"*" 1 < (-l)^(n + jfim+n-l) = a 2 

Thus the sequence {aj} is log-concave and so is {x m } due to (|3.ip and the afore- 
mentioned property of the binomial convolution. The sequence {y m } is log-concave 
because 

1 1 

ym-iym+i — JZ Z TTFTo i TTT 

(2n — 1 — to + lj! (2n — 1 — to — 1)! 

2n — 1 — to 1 1 9 



2n - 1 - m + 1 (2n - 1 - m)! (2n - 1 - to) 



□ 



Now we prove that the coefficients /3 ni fc are positive. The following two lemmas 
are useful. 

Lemma 3.2. If a positive sequence {klck} k>0 is log-concave, then 
(i) {(k + l)ck+i/ck} is non-increasing; 
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(ii) {cfc} is log-concave; 

(iii) the terms cu satisfy 

(k + to\ 

(3.2) c k c m > I JcoCfc +m (0 < to < fc + 1) . 

Proof. The statements (i) and (ii) are obvious. To prove (iii) we apply a method 
used in [1] . Specifically, by (i) we have for < p < k 

Cp+i k + p + 1 Ck+p+i 
c p ~ p + 1 c fc 

Apply the last inequality for p = 0, 1, 2, to with m < fc+1, and form the products 
of all left-hand and right-hand sides. As a result, after the cancellation we obtain 

c m > k + I k + 2 k + m c k+rn 



Co 12 to Ck 

which is equivalent to (13.21) . □ 

Lemma 3.3. If the coefficients /3 n> fc are positive, then for fixed n > 3 they satisfy 

,o q \ (fc + l)ftn,fc+l 

(3.3) ! < n - 1 . 

Pn.fe 

Proof. By Theorem 13.11 and under the assumption of lemma, for fixed n > 3 the 
sequence {fc!/?ri,fc}fc=o meet the conditions of Lemma 13.21 Applying the inequality 
(|3.2[) with to = 1 to this sequence gives (k + l)f3 n ,k+i/fln,k < Pn,i/Pn,Q- Then the 
lemma follows as due to 



— — = , = 3n — n 1 H — — 1 < 3n — 2n — 1 = n — 1 . 

□ 

Theorem 3.4. The coefficients fi n ,k ore positive. 

Proof. We prove the statement by induction on n. It is true for n < 5 (see §1). As- 
sume that for some fixed n all the members of the sequence {/3n,fc}^Zg are positive. 
Since /3 n +i,o = (n + 1)" > and /3 n+ i,„ = n! > by (j!T2"f and (|2.3p. we only need 
to consider fc = 1, 2, n — 1. 

Substituting inequalities ^„,fc+i < (n— l) J 5„ ) / s /(fc+l) and /3 n ,fc-l > k/3 n ,k/{n- 1), 
which follow from (|3.3j) . in the recurrence (|2.4|) immediately gives the result 

fc n — 1 ( fc \ 
/9n+i,ft > (3n-fc-l)/3 n . fc +n _ fc = 2n + P n ,k > • 

71 — 1 fc+1 \ 71 — 1/ 

Thus the proof by induction is complete. □ 

Corollary 3.5. The sequence {/3n,fe}fe=o * s unimodal for n > 3. 

Proof. By Theorem 13.41 the sequence {/3n,fc}fc = g is positive, therefore by Theorem 
13.11 and Lemma l3.2f ii) it is log-concave and, hence, unimodal. □ 
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4. Relation to Carlitz's numbers 

There is a relation between the coefficients f3 n ^ and numbers _B(k,j, A) intro- 
duced by Carlitz in [J]. Comparing the formulae (|2.8I) and (|5.1[) with the corre- 
sponding [H eq.(6.3)] and [H eq.(2.9)], taking into account that he uses the notation 
i?(n,m,r) = {^+;} r ,wefind 

(4.1) p nM = {-l) k B{n-l,n-\-k,n) . 

It follows that for n > 3, the sequence {Bin — 1, k, n)}2Z is log-concave together 
with {/3 n ,k}k=a- 

Using the property [4j eq.(2.7)] that J2j=o A) = (2k — 1)!!, we can compute 

p n {w) at the singular point where W = — 1 (cf. (jl.ll) ). Thus, substituting w = — 1 
in (I2.ip gives p„(— 1) = (— l) rl_1 (2n — 3)!!. Thus w = — 1 is not a zero of p n (w). 

We also note that the numbers B(K,j,X) are polynomials of A and satisfy a 
three-term recurrence [4j eq. (2.4)] 

(4.2) B(n,j, A) = ( /t + j-A) J B( K -l,i,A) + (K-i + A) J B( K -l,i-l,A) 

with B(k,0, A) = (1 — A) fe , B(0,j,X) — 5j, . This gives one more way to com- 
pute the coefficients Pn,k, specifically, for given n and k we find a polynomial 
B(n — l,n — 1 — fc, A) using (|4.2p and then set A = n to use (|4.1[) . 



5. Concluding remarks 

It has been established that the coefficients of the polynomials (— l) n ^ 1 p n (w) are 
positive, unimodal and log-concave. These properties imply an important prop- 
erty of W. In particular, it follows from formula (jl.ll) and Theorem 13.41 that 
(-l^-^dW/ete^ 71 - 1 ) > for n > 1. Since W(x) is positive for all positive x [5], 
this means that the derivative W is completely monotonic and W itself is a Bern- 
stein function [2]. 

Some additional identities can be obtained from the results above. For example, 
computing /?„.„_ i by (|2 .8[) and comparing with (|2.3I) gives 



ro=0 



2n - 1 \ (2n - 1 + m 
n — in — 1/1 n + m 



(n - 1)! 



A relation between { 2n ~+m m } n and 5;~7(n) 

can be obtained from (|2.8I) and 
but this is a special case of [H eq. (7.5)]. It is finally interesting to note that 
and (|2.9I) can be inverted. Indeed, in these formulae for fixed n, the sequence 
(— l) k Pn,k is a convolution of two sequences and the corresponding relation between 
their generating functions is G(w) = (1 — w) 2n ~ 1 F(w). Since -F(w) = G(w)(l — 



w 



(2n-i) _ q( w } £ fe>Q Cl~ 2 X k )w k , the inverse of, for example, (J2 



(2n — 1 + ml /2n-2 + m-k 

(5-1) ^ =^(-l) fc /3„ ifc 



n + m ^— ' 'V 2n — 2 

fe=o 
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